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There is an extensive history of form factor calculations on the lattice, primarily with ground states 
for both initial and final states. However, there have never been any radially excited transition 
form factor calculations. Furthermore, the lattice faces difficulty in extracting signal from noise 
at large transfer momenta (Q^). These measurements could give important theoretical input to 
experiments, such as those of JLab's 12 GeV program and studies of deformation of the nucleon. 
In this work, I will present a simple technique to resolve both of these difficulties and present 
results from anisotropic configurations showing improved signals for excited-state quantities. It 
should also be possible to apply this technique to isotropic lattices for calculating large-Q^ form 
factors. 
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1. Introduction 

Higher-2^ data will help us to understand hadrons and challenge QCD-based models. At 
Jefferson Lab, experiments are looking for whether (and if so, where) the Sachs electric form 
factor of the proton {G'^) crosses zero at large Q^. The future 12 GeV upgrade will provide data at 
even larger values of Q^. We in the lattice community would like to make predictions using QCD 
to compare with these new data sets. 

The form factors of excited states are another interesting subfield being explored. Experiments 
at Jefferson Lab (CLAS), MIT-Bates, LEGS, Mainz, Bonn, GRAAL, and Spring-8 measured he- 
licity amplitudes of the nucleon-Roper transition, and many models fail to predict its value. We 
are particularly interested in resolving the mysterious nature of the Roper, which is still not fully 
understood. Many believe that the Roper is simply the first radially excited state of the nucleon; 
we can verify this by providing the first-principles theoretical inputs to compare with experimental 
data. 

Lattice QCD has successfully calculated many hadronic form factors and provided first-principles 
theoretical inputs for exploring the structure of hadrons. However, not much progress have been 
made in transition form factors or at large Q^. The typical range in lattice calculations of 
nucleon form factors is < 3.0 GeV^; when one attempts higher-Q^ calculations, they suffer from 
poor signal-to-noise ratios [[I]]. In the case of transition form factors to radially excited states, the 
difficulty arises through the rapid exponential fall-off in Euclidean time of the excited state. By 
using anisotropic lattices, excited-state signals may be improved dramatically compared with con- 
ventional isotropic lattices. An exploratory study can be found in Ref. [^. In this work, we 
demonstrate that we can address these two lattice QCD challenges on the anisotropic lattice for the 
nucleon; the same techniques can be applied to other hadrons. 

2. Methodology and Setup 

We use quenched 16^ x 64 lattices with anisotropy ^ = 3 (that is, = 3at), using Wilson 
gauge action with j3 = 6.1 and stout-link smeared[||] Sheikholeslami-Wohlert (SW) fermions[||] 
with smearing parameters {p,«p} = {0.22,2}. The parameter v is nonperturbatively tuned using 
the meson dispersion relation, and the clover coefficients are set to their tadpole-improved values. 
The inverse spatial lattice spacing is about 2 GeV, as determined by the static-quark potential, and 
the simulated pion masses are about 480, 720 and 1 100 MeV. In total, we use 200 configurations 
at each pion mass. 

The interpolating fields we use for the nucleon have the general form for most octet baryons, 

X'{x) = e"''[qf{x)CY5qUmi^), (2-1) 

where C is the charge conjugation matrix, and qi and q2 are any of the quarks {u,d,s}. In the case 
of the proton, we want q\=u and 172 = d. 

At each pion mass, we use three Gaussian smearing parameters: a G {0.5,2.5,4.5}; the largest 
of these has excellent overlap with just the ground state, while the other two include substantial 
higher-state contributions. For both two-point and three-point correlators, we calculate all 9 possi- 
ble source-sink smearing combinations. We apply the variational method[Q] to extract the principal 
correlators corresponding to pure energy eigenstates from our matrix of correlators. This gives 3 
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Figure 1: Nucleon and Pn dispersions. The y-axis is in units of a,'^, while the x-axis is in units of jf—aj'^. 
Both figures show slopes consistent with 1 . 

principal correlators; we discard the 2"''-excited correlator which is likely contaminated by unre- 
solved higher states. We extract the energies of the ground and first-excited nucleon by applying 
simple exponential fits individually to the lowest two principal correlators. After this, we also 
check whether the original two-point correlators are recovered by appropriate combinations of the 
exponentials. Figure [l] shows the the dispersions of the nucleon and its first radially excited state 
at ruji ^ 720 MeV; they both have slopes consistent with the expected continuum value, 1 . Similar 
results also found for the other two pion masses. The excited state has a better signal compared 
with those obtained from conventional isotropic lattices. 

To calculate three-point functions, we choose source and sink locations to be at timeslices 15 
and 48 respectively. (This gives a source-sink separation about 1.1 fm.) We choose the final state 
to always have zero momentum in this proceeding, while the momentum of the initial state varies. 

The three-point function has the form 

n n' 

X Y.TapU„,{pf,s')p{N„'{pf,s')\V^\N„ipi,s))Un{pi,s)a, (2.2) 
s,s' 

where fn„i{pf, Pi, E'^,En,t,ti,tf) contains kinematic and overlap factors and time-dependent terms, 
n and n' index different energy states, the projection used is T = ^(1 + 74)(1 +27573), and Zy is 
the vector-current renormalization constant (which is set to its nonperturbative value). The vector 
current in Eq. 2.2 is 0{a) on-shell improved with the improved coefficient set to its tree-level value. 
Note that when calculating the three-point Green function in full QCD, there are two possible 
contraction topologies: "connected" and "disconnected", when the vector current vertex appears 
on a vacuum bubble. In this work, only "connected" quantities are included. 

The ratio method will not extract matrix elements beyond the ground state, but using the 
Z's and £"s derived from our analysis of the two-point functions, we can extract excited matrix 
elements from the three-point functions by fitting to the form given in Eq. 2.2. We want to keep 
terms in Eq. 2.2 having n and n' running from the ground to the first-excited state; thus, there 
would be four matrix elements in the minimal expansion, which will be free parameters in our fit. 
We increase the number of three-point correlators, first using just the diagonal correlators where 
the smearing is the same at the source and sink, and then using all 9 smearing combinations. 
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Figure 2: d quark contribution to the nucleon Dirac form factor Ff, obtained from the ratio approach (stars) 
and the fitting method (circles) 



The nucleon-nucleon matrix elements are verified against the ratio method, and we check that the 
transition matrix elements are consistent between different sets. In the rest of this work, we show 
the results from fuU 9-correlator simultaneous fits. 

The vector-current matrix elements, (A'^ | | A'^) (with n = n' = 0) and (A'^ | | A i ) (with n = 0, 
n' = 1) (or its conjugate) are fitted according to Eq. 2.2. The fit ranges are carefully chosen to give 
consistent results. The form factors and the vector-current matrix element for any nucleon states, 
A^i and N2, are connected through 



' unAp)^'''"', (2-3) 



where the equation of motion is used to simphfy —q^y^ = M^^ —M^^. By using singular-value 
decomposition (SVD), we solve for the form factors Fi^2- 

3. Numerical Results 

The ground-state nucleon form factors, F^^ and F^^, have been well studied, and we can use 
the conventional analysis method to crosscheck our fitting proposal. We use the largest Gaussian 
smearing (a = 4.5) three-point data, which has good overlap with ground states; this makes it ideal 
for use in a ratio approach|^ to get nucleon form factors. We obtain nucleon-nucleon form factors 
Fi and F2 consistent with those derived from the ratio method. Fig. ^ shows that both methods 
have consistent J-quark contributions to the Dirac form factor at pion mass 720 MeV. Note that 
we not only get consistent numbers in Fi , but at large momenta, our fitting approach dramatically 
improves the signal. This is because when one tries to improve the ground-state signal (normally 
done by examining the hadron effective-mass time dependence at zero momentum) with a choice 
of smearing operator, one wipes out not only the excited states but also the higher-momentum 
states. Therefore, when one tries to project onto these higher momenta, there is no doubt that the 
signal-to-noise ratio will worsen. Since we are explicitly considering excited states in our analysis, 
our sources which have good overlap with higher-momenta will be only lightly contaminated by 
excited-state signals. 

Figure || shows preliminary results with isovector Dirac and Pauli form factors, Fi 2 with three 
values of light-quark mass, which correspond to roughly 480 (pink triangles), 720 (purple squares) 
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Figure 3: Isovector Dirac (left) and Pauli (right) nucleon form factors, Fi and F2, with pion masses of 480 
(pink), 720 (purple) and 1 100 (grey) MeV. The dashed lines indicate the experimental values. 



and 1 100 (downward triangles) MeV pions. We see very mild pion-mass dependence over the 
entire range for the Dirac for factor and in the large- 2^ region for the Pauli form factor. Al- 
though F2 shows a small tendency toward the experimental curve, these form factors still deviate 
significantly from experiment. This is somewhat expected, since our pion masses are much heavier 
than the physical value. Even in unquenched calculations, the nucleon form factors do not agree 
with experimental values with pion masses as low as 300 MeV. (See reviews in Ref. [^] or these 
proceedings|||].) We may expect to see the lattice data approaching experiment as the pion mass 
used is decreased in future calculations. Further study below 300 MeV sea pion mass wiU be 
critical to understanding this behavior. 

Note that we get very clean signal in the range 3 GeV^ < < 5 GeV^. When we examine 
the electric Sachs form factor 



^2 



GE{q') =Fi{q')-^F2{q'), (3.1) 

we see that does not cross zero within our region. In the near future, we will continue 
analyzing data with non-zero final-state momentum to look for even larger form factors. 

We can obtain the P\i-N transition form factors from processes that are related to the Roper de- 
cay Pu ^ yN i {p\V^\Pn)) and the one that is related to photoproduction, fN ^ Ai ( (Ai | p))- 
We show preliminary results for transition form factors F]*2 at three pion masses in Figure ^ and 
Figure |5|; the experimental values are reconstructed from experimental helicity amplitudes using 
single-pion or two-pion channel analyses[^, [T^, |n|, |l2|] and PDG[ 13]. Note that since our nucleon 



and Roper masses are much higher than the physical ones, we are in the time-like region when 
we use the matrix element (Ai lA'^) to construct the transition form factors. As we decrease the 
pion mass, we will enter the space-like region, and this matrix element will be helpful in giving us 
different Q^. 

A Pi 1 at rest on the lattice has a decay into a pion and a nucleon in a f -wave and a decay into 
two pions and a nucleon in an S-wave via ^n^r. Since we do not wish to consider the complicated 
case that occurs when multi-particle states may be present, we must avoid kinematical situations in 
which decays can occur. In our simulation, these are suppressed by the high quark mass and dis- 
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Figure 4: Proton-Pn transition form factors: 
(pink), 720 (purple) and 1100 (grey) MeV. The brown points are the experimental results from CLAS and 
the orange are numbers from the PDG. Two matrix elements are used to extract the form factors: (p | | A i ) 
(diamonds) and {Pn \ V^\p) (squares). 



cretization of momentum. However, the lowest-momentum Pi i can decay into a lowest-momentum 
pion and a nucleon at rest. Since we cannot untangle the Pi i from the two-particle state here, the 
leftmost points from (p | | Pn) elements should be disregarded. 

In the time-like regions, all four transition form factors have very mild quark-mass depen- 
dence. Even in the space-like regions, the dependence only become evident at 2^ > 2 GeV^. We 
obtain form factors with momentum transfer up to 6 GeV^; unfortunately, our lightest pion-mass 
calculation suffers from large statistical error; we need more statistics to give a clean signal at 
large Q^. In the large-Q^ region, our calculations show that these form factors have about the same 
magnitude with the single-pion or two-pion analyses of the experimental data. However, in the 
low-Q^ region, both the CLAS analysis and PDG value (at zero momentum transfer) have opposite 
sign from ours for both proton and neutron transition form factors. Pion cloud effects could be 
significantly larger in the low region; therefore, removal of the quenched approximation would 
be very helpful for distinguishing whether there is a true discrepancy in the sign. 

4. Conclusion 

In this work, we try to solve two form-factor challenges in lattice QCD calculations: large-Q^ 
form factors and excited-ground state transition form factors. We demonstrate the method using 
a baryonic system, which typically has worse signal-to-noise ratio than mesons. The method can 
easily be applied to other hadronic systems. 

Using multiple operators to determine the ground and excited states, we are able to extract 
good signals for both the ground-state form factors at large and the excited-ground transition 
form factors. In this work, we construct a basis from operations with various smearing parameters. 
We see significantly improved signal-to-noise ratios for large-Q^ momentum N-N form factors 
compared with previous studies. We find no zero-crossing for up to 5.5 GeV^. We calculate 
N — Pii transition form factors up to around 6 GeV^ and find small quark-mass dependence. We 
would like to move on to dynamical calculations to re-examine the low-2^ behavior. Future gener- 
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Figure 5: Neutron-Pn transition form factors: F" ^" (left) and F2 (right) with pion masses of 480 
(pink), 720 (purple) and 1 100 (grey) MeV. The symbols are the same as in Figure 4. 



alizations to operators constructed in irreducible representations of the cubic group will probably 
allow us to analyze form factors from more than one radially excited state. 
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